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Abstract 

For an equiregular sub-Riemannian manifold M, Popp's volume is a smooth volume 
which is canonically associated with the sub-Riemannian structure, and it is a natural gen- 
eralization of the Riemannian one. In this paper we prove a general formula for Popp's 
volume, written in terms of a frame adapted to the sub-Riemannian distribution. As a 
hrst application of this result, we prove an explicit formula for the canonical sub-Laplacian, 
namely the one associated with Popp's volume. Finally, we discuss sub-Riemannian isome- 
tries, and we prove that they preserve Popp's volume. We also show that, under some 
hypotheses on the action of the isometry group of M, Popp's volume is essentially the 
unique volume with such a property. 

1 Introduction 

The problem to define a canonical volume on a sub-Riemannian manifold was first pointed out 
by Brockett in his seminal paper [10] . motivated by the construction of a Laplace operator on 
a 3D sub-Riemannian manifold canonically associated with the metric structure, analogous to 
the Laplace-Beltrami operator on a Riemannian manifold. Recently, Montgomery addressed 
this problem in the general case (see [131 Chapter 10]). 

Even on a Riemannian manifold, the Laplacian (defined as the divergence of the gradient) is 
a second order differential operator whose first order term depends on the choice of the volume 
on the manifold, which is required to define the divergence. Naively, in the Riemannian case, 
the choice of a canonical volume is determined by the metric, by requiring that the volume of 
a orthonormal parallelotope (i.e. whose edges are an orthonormal frame in the tangent space) 
is 1. 

From a geometrical viewpoint, sub-Riemannian geometry is a natural generalization of Rie- 
mannian geometry under non-holonomic constraints. Formally speaking, a sub-Riemannian 
manifold is a smooth manifold M endowed with a bracket-generating distribution T> C TM, 
with k = rankD < n = dimM, and a smooth fibre- wise scalar product on T>. From this struc- 
ture, one derives a distance on M - the so-called Carnot-Caratheodory metric - as the infimum 
of the length of horizontal curves on M, i.e. the curves that are almost everywhere tangent to 
the distribution. 

Nevertheless, sub-Riemannian geometry enjoys major differences with respect to the Rie- 
mannian case. For instance, a construction analogue to the one described above for the Rie- 
mannian volume is not possible. Indeed the inner product is defined only on a subspace of the 
tangent space, and there is no canonical way to extend it on the whole tangent space. 
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Popp's volume is a generalization of the Riemannian volume in sub-Riemannian setting. It 
was first denned by Octavian Popp but introduced only in [13] (see also [3]). Such a volume 
is smooth only for an equiregular sub-Riemannian manifold, i.e. when the dimensions of the 
higher order distributions V 1 := V, V i+1 := + [V\V], for every i > 1, do not depend on the 
point (for precise definitions, see Sec. 2). 

Under the equiregularity hypothesis, the bracket-generating condition guarantees that there 
exists a minimal m £ N, called step of the structure, such that T> m = TM. 

Then, for each q £ M, it is well defined the graded vector space: 

m 

gr q (V) := V l q /V\-\ where V° q = 0. (1) 

i=i 

The vector space gr„(P), which can be endowed with a natural sub-Riemannian structure, is 
called the nilpotentization of the structure at the point q, and plays a role analogous to the 
Euclidean tangent space in Riemannian geometry. Popp's volume is defined by inducing a 
canonical inner product on gr q (T>) via the Lie brackets, and then using a non- canonical isomor- 
phism between gr (P) and T q M to define an inner product on the whole T q M. Interestingly, 
even though this construction depends on the choice of some complement to the distribution, 
the associated volume form (i.e. Popp's volume) is independent on this choice. 

It is worth to recall that on a sub-Riemannian manifold, which is a metric space, the Hauss- 
dorff volume and the spherical Hausdorff volume, respectively and S®, are canonically de- 
finedQ The relation between Popp's volume and has been studied in [2], where the authors 
show how the Radon-Nikodym derivative is related with the nilpotentization of the structure. 
In particular they prove that the Radon-Nikodym derivative could also be non smooth (see also 
[6j[8]). Remember that the Hausdorff and spherical Hausdorff volumes are both proportional to 
the Riemannian one on a Riemannian manifold. The relation between Hausdorff measures for 
curves and different notions of length in sub-Riemannian geometry is also investigated in |11| . 

On a contact sub-Riemannian manifold, Popp's volume coincides with the Riemannian vol- 
ume obtained by "promoting" the Reeb vector field to an orthonormal complement to the 
distribution. In the general case, unfortunately, the definition is more involved. To the authors' 
best knowledge, explicit formulas for Popp's volume appeared, for some specific cases, only in 

mum- 

The goal of this paper is to prove a general formula for Popp's volume, in terms of any 
adapted frame of the tangent bundle. In order to present the main results here, we briefly 
introduce some concepts which we will elaborate in details in the subsequent sections. Thus, 
we say that a local frame X\, . . . , X n is adapted if X%, . . . , X^. is a local frame for P l , where 
ki := dimP*, and X\, . . . ,Xj~ are orthonormal. Even though it is not needed right now, it is 
useful to define the functions c\j G C°°{M) by 

n 

[X i ,X j ] = J24 j Xi. (2) 
1=1 

With a standard abuse of notation we call them structure constants. For j = 2, . . . , m we define 
the adapted structure constants b[ i ^ € C°°(M) as follows: 

[X il ,[X i2 ,...,[X ij _ 1 ,X ij }]}= /'!,,,..., .V, modi*'- 1 , (3) 

l=kj-!+l 

1 Recall that the Hausdorff dimension of a sub-Riemannian manifold M is given by the formula Q = X^fci * n «i 
where rij := dim2?*/X>* -1 . In particular the Hausdorff dimension is always bigger than the topological dimension. 
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where 1 < < k. These are a generalization of the c-.-, with an important difference: 

J LJ 

the structure constants of Eq. ([2]) are obtained by considering the Lie bracket of all the fields of 
the local frame, namely 1 <i,j,l < n. On the other hand, the adapted structure constants of 
Eq. are obtained by taking the iterated Lie brackets of the first k elements of the adapted 
frame only (i.e. the local orthonormal frame for D), and considering the appropriate equivalence 
class. For j = 2, the adapted structure constants can be directly compared to the standard 
ones. Namely b\j = c\j when both are defined, that is for 1 < i, j < k, I > k + 1. 
Then, we define the kj — kj-\ dimensional square matrix Bj as follows: 

k 

[Bjf= £ l>l,..,l>,, : ,..,,- J 1 »<• (4) 

with the understanding that B\ is the kx k identity matrix. It turns out that each Bj is positive 
definite. 

Theorem 1. Let X\,... ,X n be a local adapted frame, and let u\ . . . ,u n be the dual frame. 
Then Popp's volume V satisfies 

V = . 1 v 1 A ... A v n , (5) 

where Bj is defined by Q in terms of the adapted structure constants ([3]) . 

To clarify the geometric meaning of Eq. ([5]), let us consider more closely the case m = 2. 
If V is a step 2 distribution, we can build a local adapted frame {X\, . . . , Xk, Xk+i, . . . ,X n } 
by completing any local orthonormal frame {X\, . . . , X^} of the distribution to a local frame of 
the whole tangent bundle. Even though it may not be evident, it turns out that B^ (q) is the 
Gram matrix of the vectors X^+i, . . . , X n , seen as elements of T q M/D q . The latter has a natural 
structure of inner product space, induced by the surjective linear map [ , ] : T> g (8) T> q — )■ T q M/V q 
(see Lemma [T3]) . Therefore, the function appearing at the beginning of Eq. ([5]) is the volume 
of the parallelotope whose edges are X\ , . . . , X n , seen as elements of the orthogonal direct sum 
gv q {V)=V q ®T q M/V q . 

With a volume form at disposal, one can naturally define the associated divergence opera- 
tor, which acts on vector fields. Moreover, the sub-Riemannian structure allows to define the 
horizontal gradient of a smooth function. Then, we define a canonical sub-Laplace operator as 
A := div o V, which generalizes the Laplace-Beltrami operator. This is a second order differ- 
ential operator, which has been studied in [31 [5]. As a corollary to Theorem [TJ we obtain a 
formula for the sub-Laplacian A in terms of any local adapted frame. 

Corollary 2. Let X\, . . . ,X n be a local adapted frame. Let A be the canonical sub-Laplacian. 
Then 

k I 1 m n \ 

A = E X? - 2 E ^(BfX^Bj)) + (6) 
i=l \ j=l 1=1 J 

where c\a are the structure constants ([2]), and Bj is defined by Q in terms of the adapted 
structure constants ([3]). 

Remark 1. If M is a Carnot group (i.e. a connected, simply connected nilpotent group, whose 
Lie algebra is graded, and whose sub-Riemannian structure is left invariant) the Bj are constant. 
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Moreover, Vi Yl?=i c \i = 0, as a consequence of the graded structure. Then, in this case, the 
sub-Laplacian is a simple "sum of squares" A = ^2i=i-^f- This is a manifestation of the fact 
that Carnot groups are to sub-Riemannian geometry as Euclidean spaces are to Riemannian 
geometry. Indeed, on M. n , the Laplace-Beltrami operator is a simple sum of squares. 

More in general, in [3], the authors prove that for left-invariant structures on unimodular 
Lie groups the sub-Laplacian is a sum of squares. 

In the last part of the paper we discuss the conditions under which a local isometry preserves 
Popp's volume. In the Riemannian setting, an isometry is a diffeomorphism such that its 
differential is an isometry for the Riemannian metric. The concept is easily generalized to the 
sub-Riemannian case. 

Definition 3. A (local) diffeomorphism (j) : M — > M is a (local) isometry if its differential 
0* : TM —> TM preserves the sub-Riemannian structure (T>, (■,■)), namely 

i) fa{V q )=V m for all q G M, 

ii) (faX, faY) Hq) = {X, Y) q for all q G M, X, Y G V q . 

Remark 2. Condition i), which is trivial in the Riemannian case, is necessary to define isometries 
in the sub-Riemannian case. Actually, it also implies that all the higher order distributions are 
preserved by fa, i.e. fa(V q ) = T^ % ^ q y for 1 < i < m. 

Definition 4. Let M be a manifold equipped with a volume form [i G £l n (M). We say that a 
(local) diffeomorphism <p '■ M — > M is a (local) volume preserving transformation if fafj, = \i. 

In the Riemannian case, local isometries are also volume preserving transformations for the 
Riemannian volume. Then, it is natural to ask whether this is true also in the sub-Riemannian 
setting, for some choice of the volume. The next proposition states that the answer is positive 
if we choose Popp's volume. 

Proposition 5. Sub-Riemannian (local) isometries are volume preserving transformations for 
Popp's volume. 

Proposition [5] may be false for volumes different than Popp's one. We have the following. 

Proposition 6. Let Iso(M) be the group of isometries of the sub-Riemannian manifold M . If 
Iso(M) acts transitively on M, then Popp's volume is the unique volume (up to multiplication 
by scalar constant) such that Proposition^ holds true. 

Definition 7. Let M be a Lie group. A sub-Riemannian structure (M,T>, (-, •)) is left invariant 
if Vg G M, the left action L g : M — > M is an isometry. 

As a trivial consequence of Proposition [S] we recover a well-known result (see again |13j). 

Corollary 8. Let (M, T>, (•,•)) be a left-invariant sub-Riemannian structure. Then Popp's 
volume is left invariant, i.e. L* g V = V for every g G M. 

Propositions [5l [6] and Corollary [8] should shed some light about which is the "most natural" 
volume for sub-Riemannian manifold. 
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2 Sub-Riemannian geometry 



We recall some basic definitions in sub-Riemannian geometry. For a more detailed introduction, 

see nnang. 

Definition 9. A sub-Riemannian manifold is a triple (M,T>, (•, •)), where 
(i) M is a connected orientable smooth manifold of dimension n > 3; 
(ii) T> C TM is a smooth distribution of constant rank k < n; 

(Hi) (•, -) g is an inner product on the fibres T> q , smooth as a function of q. 

Let T(T>) C Vec(M) be the C°°(M)-module of the smooth sections of T>. Throughout 
this paper we assume that the sub-Riemannian manifold M satisfies the bracket- generating 
condition, i.e. 

span{[Xi, [X 2 ,..., [X.^XjMq) | X t E T(V), j G N} = T,M, Vg£M. (7) 

In other words, the iterated Lie brackets of smooth sections of T> span the whole tangent bundle 
TM. Condition (|7|) is also called Hormander condition, and bracket-generating distribution are 
also referred to as completely nonholonomic distributions. 

An absolutely continuous curve 7 : [0, T] — )■ M is said to be horizontal (or admissible) if 

j(t) G P 7(t ) for a.e. f E [0,7]. 

Given an horizontal curve 7 : [0, T] —> M, the length 0/7 is 

%)= TitWI d*. (8) 

JO 

The distance induced by the sub-Riemannian structure on M is the function 

d(qo,qi) = inf{^(7) I 7(0) = qo,l(T) = gi,7 horizontal}. (9) 

The connectedness of M and the bracket-generating condition guarantee the finiteness and the 
continuity of the sub-Riemannian distance with respect to the topology of M (Chow-Rashevsky 
Theorem, see, for instance, [1]). The function d(-,-) is called the Carnot-Caratheodory distance 
and gives to M the structure of metric space (see [3 Q2] ). 

Locally (i.e. on an open set U C M), there always exists a set of k smooth vector fields 
X\, . . . , Xk such that, \/q G U, it is an orthonormal basis of V q . The set {X\, . . . , X^} is called 
a local orthonormal frame for the sub-Riemannian structure. 

Definition 10. Let P be a distribution. Its flag at q G M is the sequence of vector spaces 
^CDjc^C.C T q M defined by 

V\ := {0}, Z>J := V q , V\ +1 := V\ + [V\V] q , 

where, with a standard abuse of notation, we understand that [D % ,T>] q is the vector space 
generated by the iterated Lie brackets, up to length i, of local sections of the distribution, 
evaluated at q. 
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Even though the rank of T> is constant, the dimensions of the subspaces of the flag, i.e. 
the numbers ki(q) := dim(P*) may depend on the point. Observe that the bracket-generating 
condition can be rewritten as follows 

Vg G M 3 minimal m{q) G N such that k m (q) = dim T q M. 

The number m(q) is called the step of the distribution at the point q. The vector G(q) '■= 
(ki(q), &2 (</),... , k m (q)) is called the growth vector of the distribution at q. 

Definition 11. A distribution D is equiregular if the growth vector is constant, i.e. for each 
i = 1,2, ... ,m, h(q) = dim(P*) does not depend on q G M. In this case the subspaces V q are 
fibres of the higher order distributions T> % C TM. 

For equiregular distributions we will simply talk about growth vector and step of the distri- 
bution, without any reference to the point q. 

Finally, we introduce the nilpotentization of the distribution at the point q, which is funda- 
mental for the definition of Popp's volume. 

Definition 12. Let T> be an equiregular distribution of step m. The nilpotentization of V at 
the point q G M is the graded vector space 

gr g (v) = v q @ v 2 jv q e . . . e v™/v™~ x . 

The vector space gx q (T>) can be endowed with a Lie algebra structure, which respects the 
grading. Then, there is a unique connected, simply connected group, Gi q {V), such that its 
Lie algebra is gr q (T>). The global, left-invariant vector fields obtained by the group action on 
any orthonormal basis of T> q C gr q (T>) defines a sub-Riemannian structure on Gr 9 (D), which is 
called the nilpotent approximation of the sub-Riemannian structure at the point q. 



3 Popp's volume 

In this section we provide the definition of Popp's volume, and we prove Theorem [TJ Our 
presentation follows closely the one that can be found in |13j . The definition rests on the 
following lemmas. 

Lemma 13. Let E be an inner product space, and let ir : E — )• V be a surjective linear map. 
Then ir induces an inner product on V such that the length of v G V is 

\\v\\v = min{ 1 1 e 1 1 £; s.t. 7r(e) = v} . (10) 

Proof. It is easy to check that Eq. (|10|) defines a norm on V. Moreover, since || • \\e is induced 
by an inner product, i.e. it satisfies the parallelogram identity, it follows that || • \\y satisfies the 
parallelogram identity too. Notice that this is equivalent to consider the inner product on V 
defined by the linear isomorphism tt : (ker-7r)^ — > V. Indeed the length of v G V is the length 
of the shortest element e G 7r _1 (t>). □ 

Lemma 14. Let E be a vector space of dimension n with a flag of linear subspaces {0} = F° C 
F 1 C F 2 C ... C F m = E. Let gr(F) = F 1 F 2 /F 1 ... F m /F m - 1 be the associated graded 
vector space. Then there is a canonical isomorphism 9 : A n E — > A n gr(F). 
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Proof. We only give a sketch of the proof. For < i < m, let ki := dimF 1 . Let X\, . . . ,X n be a 
adapted basis for E, i.e. X\, . . . , X^ is a basis for i 7 *. We define the linear map 9 : E — > gr(-F) 
which, for < j < m — 1, takes -Xj^+i, . . . ,X^. +1 to the corresponding equivalence class in 
F^ +1 j ' FK This map is indeed a non-canonical isomorphism, which depends on the choice of 
the adapted basis. In turn, 9 induces a map 9 : A n E — > A n gr(F), which sends X± A . . . A X n 
to 9{X\) A ... A 0(X n ). The proof that 9 does not depend on the choice of the adapted basis is 
"dual" to [13 Lemma 10.4]. □ 

The idea behind Popp's volume is to define an inner product on each T> q /T> % ~ 1 which, in 
turn, induces an inner product on the orthogonal direct sum gr q (P). The latter has a natural 
volume form, which is the canonical volume of an inner product space obtained by wedging 
the elements an orthonormal dual basis. Then, we employ Lemma [TH to define an element of 
(f\ n T q M)* ~ A n T*M, which is Popp's volume form computed at q. 

Fix q G M. Then, let v, w G T> q , and let V, W be any horizontal extensions of v, w. Namely, 
V, W E T(V) and V(q) = v, W(q) = w. The linear map vr : V q <g> V q -> V%/V q 

tt(v®w) :=[V,W] q modP 9 , (11) 

is well defined, and does not depend on the choice the horizontal extensions. Indeed let V and 
W be two different horizontal extensions of v and w respectively. Then, in terms of a local 
frame X\ , . . . , X^ of T> 

k k 

V = V + ^f l X l , W = W + Y J 9iX l , (12) 
i=l i=l 

where, for 1 < i < k, fi,gi G C°°(M) and fi(q) = gi(q) = 0. Therefore 

k k 

[V, W] = [V, W] + Y, - W Ui)) X i + E foil*,*}] ■ (13) 



Thus, evaluating at q, [V,W] q = [V,W] q mod V q , as claimed. Similarly, let 1 < i < m. The 



linear maps 7Tj : &V q -^V q /V q 1 



7Ti(vi (8) • • • (8 t><) = [Fi, [V^, . . . , [Vi- U Vi\]] q mod V\~ x , (14) 

are well defined and do not depend on the choice of the horizontal extensions V\ , . . . , Vi of 
vi, ...,«». 

By the bracket-generating condition, 7Tj are surjective and, by Lemma [T3| they induce an 
inner product space structure on Therefore, the nilpotentization of the distribution 

at q, namely 

g r q (v) = v q ® v 2 q /v q e . . . e v^/v™- 1 , (15) 

is an inner product space, as the orthogonal direct sum of a finite number of inner product 
spaces. As such, it is endowed with a canonical volume (defined up to a sign) fi q G A n gr (? (D)*, 
which is the volume form obtained by wedging the elements of an orthonormal dual basis. 

Finally, Popp's volume (computed at the point q) is obtained by transporting the volume 
of gr q (D) to T q M through the map 9 q : f\ n T q M —> A n gr q (D) defined in Lemma [T4J Namely 

v q = e*M = N oe q , (i6) 
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where 6* denotes the dual map and we employ the canonical identification (A n T q M)* ~ A n T*M. 
Eq. ()16p is defined only in the domain of the chosen local frame. Since M is orientable, with a 
standard argument, these n-forms can be glued together to obtain Popp's volume V G Q n (M). 
The smoothness of V follows directly from Theorem [TJ 

Remark 3. The definition of Popp's volume can be restated as follows. Let (M, D) be an oriented 
sub-Riemannian manifold. Popp's volume is the unique volume V such that, for all q G M, the 
following diagram is commutative: 

(M,V) (A n T g M)* 

gx q (V) > (A n gr q (V))* 

where /i associates the inner product space gr q (V) with its canonical volume fi q , and 6* is the 
dual of the map defined in Lemma HH 



3.1 Proof of Theorem Q] 

We are now ready to prove Theorem[TJ For convenience, we first prove it for a distribution of step 
m = 2. Then, we discuss the general case. In the following subsections, everything is understood 
to be computed at a fixed point q G M. Namely, by gr(P) we mean the nilpotentization of T> 
at the point q, and by T> 1 we mean the fibre T> l q of the appropriate higher order distribution. 

3.1.1 Step 2 distribution 

If P is a step 2 distribution, then T> 2 = TM. The growth vector is Q = (k, n). We choose n — k 
independent vector fields {Yi}f =k+1 such that X\, . . . , X^, ifc+i, . . . , Y n is a local adapted frame 
for TM. Then 

n 

[x i ,x j ]= h % Y l mod ^- ( 17 ) 

l=k+l 

For each I = k + 1, . . . , n, we can think to b\j as the components of an Euclidean vector in IR fc2 , 
which we denote by the symbol b l . According to the general construction of Popp's volume, we 
need first to compute the inner product on the orthogonal direct sum gr(P) = T> © T> 2 /T>. By 
Lemma [T31 the norm on T> 2 /T> is induced by the linear map tt : ® 2 T> — > T> 2 /T> 

ir(Xi (8) Xj) = [Xi,Xj] mod V . (18) 

The vector space (3 2 T> inherits an inner product from the one on V, namely \/X, Y, Z,W G V, 
(X®y, Z®W) = (X, Z)(Y, W). 7r is surjective, then we identify the range V 2 /V with ker tt 1 - C 
(& 2 T>, and define an inner product on T> 2 /T> by this identification. In order to compute explicitly 
the norm on T) 2 I'D (and then, by polarization, the inner product), let Y G T> 2 /T>. Then 

\\Y\\ V2/V = min{||Z||^ s.t. n{Z) = Y} . (19) 

Let Y = J2?=k+i and Z = ^ - =1 aijXi®Xj G ® 2 V. We can think to a>ij as the components 
of a vector a G M fc2 . Then, Eq. (|19|) writes 

II^IId 2 /© = min{|a| s.t. a ■ b l = c\ I = k + 1, . . . , n} , (20) 
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where \a\ is the Euclidean norm of a, and the dot denotes the Euclidean inner product. Indeed, 
ll^llx> 2 I'D is the Euclidean distance of the origin from the affine subspace of M k defined by the 
equations a ■ b l = c l for I = k + 1, . . . , n. In order to find an explicit expression for ||y||^ja m m 
terms of the b 1 , we employ the Lagrange multipliers technique. Then, we look for extremals of 

n 

L(a, b k+1 , ...,&«, X k+1 , . . . , A n ) = \a\ 2 - 2 £ X^a ■ b l - c l ) . (21) 

l=k+l 

We obtain the following system 

n 

\ r b l -a = 0, 



=k+l 
n 

Xib 1 -b r = c r , r = k + l,...,n. 



(22) 



Let us define the n — k square matrix B, with components B hl = b h ■ b l . B is a Gram matrix, 
which is positive definite iff the b l are n — k linearly independent vectors. These vectors are 
exactly the rows of the representative matrix of the linear map tt : ® 2 V -^V 2 /V, which has 
rank n — k. Therefore B is symmetric and positive definite, hence invertible. It is now easy to 
write the solution of system (|22p by employing the matrix B" 1 , which has components B7, , 
Indeed a straightforward computation leads to 

||c s y s |||,2/x) = ^B^c 1 . (23) 

By polarization, the inner product on T> 2 /T> is defined, in the basis Yi, by 

{Y h Y h ) v * /v = Bi h l . (24) 

Observe that B~ l is the Gram matrix of the vectors . . . ,Y n seen as elements of T> 2 /V. 

Then, by the definition of Popp's volume, if u 1 , . . . , v k , jU fc+1 , . . . , [i n is the dual basis associated 
with Xi, . . . , Xk, Yk+i, . . . ,Y n , the following formula holds true 

V = . 1 v 1 A • • • A v k A fi k+1 A • • • A /Z 1 . (25) 
vdet B 

3.1.2 General case 

In the general case, the procedure above can be carried out with no difficulty. Let X%, . . . ,X n 
be a local adapted frame for the flag V° C V C V 2 C • • • C V m . As usual ki = dim(P*). For 
j = 2, . . . ,m we define the adapted structure constants b\ i i . £ C°°(M) by 

[X h ,[X i2 ,...,[X ij _ l ,X ij }]]= mod^- 1 , (26) 

l=kj -i + l 

where 1 < i\, . . . , ij < k. Again, b\ x i . can be seen as the components of a vector b l G M fcJ . 
Recall that for each j we defined the surjective linear map irj : &T> — > T>i /T>i~ l 

tt, {X h ®X i2 ®---®Xi.) = [X H ,[X i2 ,..., [X^ , Xi .]]] mod Vi' 1 . (27) 
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Then, we compute the norm of an element of T>i /X*- 7-1 exactly as in the previous case. It is 
convenient to define, for each 1 < j < m, the kj — kj-i dimensional square matrix Bj, of 
components 

m hl = e b U,.,A^y ( 28 ) 

ii,i2,...,ij=l 

with the understanding that B\ is the k x k identity matrix. Each one of these matrices 
is symmetric and positive definite, hence invertible, due to the surjectivity of ttj. The same 
computation of the previous case, applied to each X> J ' jV^~ x shows that the matrices BJ 1 arc 
precisely the Gram matrices of the vectors X kj _ 1 +i, . . . , X kj £ P- 7 /X> 3 ' -1 , in other words 

(X kj _ 1+u X kj _ 1+h ) VJ /V] -i = Bg 1 . (29) 

Therefore, if v 1 , . . . , v n is the dual frame associated with X\, . . . , X n , Popp's volume is 

V = . 1 v 1 A . . . A u n . (30) 

3.2 Examples 

In this section we compute Popp's volume for some specific equiregular sub-Riemannian struc- 
tures. We also discuss, through an example, the non-equiregular case. 

3.2.1 Contact manifolds 

Contact manifolds are a well-known class of sub-Riemannian structures. We recall the basic 
definition first, then we compute Popp's volume in terms of a canonical operator associated 
with a contact structure. 

Definition 15. Let co G Q l (M) be a one-form on M. Let V be the n— 1 dimensional distribution 
T> := kerw. We say that uj is a contact form if <&j|d is non degenerate. In this case, D is a 
called contact distribution. A sub-Riemannian structure (M,T>, (•,•)), where 2? is a contact 
distribution, is is called contact sub-Riemannian manifold. 

Notice that the non-degeneracy assumption implies that the dimension of M is odd. Observe 
that any contact manifold satisfies the bracket-generating condition, is equiregular, has step 2, 
and its growth vector is Q = (n — 1, n). 

With any contact form oj we can associate its Reeb vector field, which is the unique vector 
field Xq satisfying the conditions uj(Xq) = 1 and dw(Xo, •) = 0. Notice that, given a local 
orthonormal frame X±, . . . , X k for the distribution, then Xi, . . . , X^, Xq is a local adapted frame, 
since Xq is transversal to T>. 

The contact form u induces a linear bundle map (i.e. a fibre-wise linear map) J : D — > T>, 
defined by (JX,Y) = dw(X,Y), VX,Y € V. Observe that the restriction J q of J to the fibres 
of V is a linear skew-symmetric operator on the inner product space (D q ,(-,-) q ). Hence its 
Hilbert-Schmidt norm \\J q \\ is well defined by the formula || J q \\ 2 = Ylij=i(Xi, JXj} 2 . 

Proposition 16. Let M be a contact sub-Riemannian manifold and J : V — > V as above. Let 
v 1 , . . . , v k , u° be the dual frame associated with the local adapted frame X\ , . . . , X k , Xq , where 
Xq is the Reeb vector field. Then 

V = -r^—V 1 f\... f\V k f\V° , (31) 

where \\J q \\ is the Hilbert-Schmidt norm of J q . 
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Proof. Let X±, . . . , X^, Xq be a local adapted frame, where Xq is the Reeb vector field associated 
with the contact form. Then, for 1 < i, j < k, the structure constants satisfy 

k 

[X h X j ] = Y,c\ j X l + c%X Q , (32) 
1=1 

k 

[X i ,X ]=J24oXi- (33) 
i=i 

By Eq. fl5}, V = v l A . . . A u k A v° where g = 1/ Y^i,j=i( c ij) 2 ■ Then the statement follows 
from the identity 

k k k k 

p\\ 2 = j^ix^jxtf = MX,*,) 2 = E ^[^ x 3]f = E(4) 2 - (34) 

i,j=i *j=i *j=i *)i=i 

Observe that, in the last equality of Eq. (|34p . we employed Cartan formula for the differential 
of a one-form, and the fact that uj(Xi) = 0. □ 

Eq. (|3"Tj) can be expressed in terms of the eigenvalues of J. See also [2J Remark 30], where 
the authors exhibit this formula for the case Q = (4,5). 

Remark 4. Let / G C°°(M) be a smooth, non- vanishing function. Then u and u ! := foj define 
the same contact distribution D. However duj' ^ fdoj and, in general, the associated Reeb 
vector field is different. On the other hand, as a consequence of the identity du'\x> = fdoo\x>, it 
follows that J' = fJ. Therefore, it is convenient to choose a "normalized" contact form, which 
is uniquely specified (up to a sign) by the condition || J q \\ 2 = 1, \/q G M. Then, in terms of the 
Reeb vector field associated with the normalized contact form, V = v 1 A . . . A u k A u°. 



3.2.2 Carnot groups of step 2 

A Carnot group G of step 2 is a left-invariant sub-Riemannian structure on a nilpotent, con- 
nected, simply connected Lie group whose Lie algebra g admits a stratification q = V\ ® V2 with 
[XijVi] = V2 and [Vi,V2] = [V^,^] = {0}. The sub-Riemannian structure is defined by left 
translation of the subspace V% , where we choose an orthonormal basis X\ , . . . , Xk ■ It is possible 
to choose a basis l^+i, . . . , Y n of Vi such that 

n 

[Xi,Xj]= Y, h % Y ^ [X i ,Y h ] = [Y h ,Y l }=0. 

h=k+l 

Using the standard exponential coordinates (i.e. the identification of the Lie group and its Lie 
algebra via the exponential map) the explicit expression for the associated left-invariant vector 
fields in R n = {(x, y) \ x G R k , y G R n - k } is 

Xi = d x . + -J2b^Xjd yh , i = l,...,k, (35) 

Y h = d Vh , h = k + l,...,n. (36) 

In [6], the authors employed the skew-symmetric matrices L h , k + 1 < h < n, of components 
[L h ]ij = b^j in order to investigate the nilpotent approximation of a step 2 sub-Riemannian 
structure. In terms of these matrices, 

B hl = {L h ,L l ), (37) 
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where (M, N) := Tr(M T N) is the Hilbert-Schmidt inner product on GL(/c, R). If the L matrices 
are orthonormal, Eq. ([5]) gives 

V = dx 1 A ... A dx k A dy k+1 A ... A dy n . (38) 

The last formula is (up to a constant factor) the definition of Popp's volume employed in [6j 
Definition 4] and [8], given in terms of a global adapted frame. 



3.2.3 Non-equiregular case 

The basic example of a bracket-generating, non-equiregular sub-Riemannian structure is the so- 
called Martinet distribution. This is the distribution on R 3 defined by the kernel of the one-form 
9 := dz — y 2 dx. A global frame for V, which we declare orthonormal, is 

X = d x + y 2 d z , Y = d y . (39) 

Let Z := d z . Then [X,Y] = -2yZ and [Y, [X,Y\] = 2Z. Observe that X,Y,Z is a global 
(adapted) frame for TM, therefore Martinet distribution is bracket-generating. However, its 
growth vector is 

[(2,2,3) ify = 0, 

and the distribution is not equiregular on the hyperplane y = 0. Nevertheless, if we restrict to 
the connected components of {y / 0}, we obtain a step 2 equiregular sub-Riemannian manifold. 
Here, Theorem [T] gives the following expression: 

V = —dx Ady Adz. (41) 

\y\ 

Eq. (|4ip shows that singularities arise precisely on the hypersurface where the equiregularity 
hypotesis fails. In [9], the authors investigate the properties of the sub-Laplacian associated 
with this volume in the Martinet structure. They show that the sub-Laplacian is essentially 
self-adjoint in each connected component of {y ^ 0}, hence the hyperplane {y = 0} acts as a 
barrier for the heat propagation. 



4 Sub-Laplacian 

In this section we define the canonical sub-Laplacian associated with a generic volume form and 
we prove Corollary [21 namely an explicit formula for the sub-Laplacian associated with Popp's 
volume. 

On a Riemannian manifold, the Laplace-Beltrami operator is defined as the divergence of 
the gradient. This definition can be easily generalized to the sub-Riemannian setting. 

Definition 17. Let / G C°°(M). The horizontal gradient of / is the unique horizontal vector 
field V/ such that 

(Vf,x) = x(f), vier(p). (42) 

It follows from the definition that, in terms of a local frame X±, . . . , X^ for T> 

k 

Vf = J2 X i(f)*i- (43) 

i=l 
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Definition 18. Let \i G O n (M) be a positive volume form, and X G Vec(M). The /i- divergence 
of X is the smooth function div^X defined by 

Lxli = div^A> . (44) 

where Cx is the Lie derivative in the direction X. 

Notice that the definition of divergence does not depend on the orientation of M, namely the 
sign of \i. The divergence measures the rate at which the volume of a region changes under the 
integral flow of a field. Indeed, for any compact C M and t sufficiently small, let e tx : O — > M 
be the flow of X G Vec(M), then 

d 
df 



fi = — I div M A> . (45) 
t=0 Je tx (n) Jn 



The next proposition is sometimes employed as an alternative definition of divergence. Let 
C^°(M) be the space of smooth functions with compact support. 

Proposition 19. For any f G Cfi°(M) and X G Vec(M) 



fdiv lx X» = - / X{f)t*. (46) 

M JM 

Proof. The proof is an easy consequence of the definition of //-divergence. □ 

The next lemma gives the relation between divergences associated with different volumes. 

Lemma 20. Let fj,, // G Q n (M) be volume forms. Let f G C°°(M), / / such that fjf = ffj,. 
Then, for any X G Vec(M) 

dhyX = div^X + X(log /) . (47) 
Proof. It follows from the Leibniz rule Cx(fn-) = {Xf)n + fCx^ = {X(\og f) + div^X)fn. □ 

When no confusion may arise, we write "div", without any reference to the volume form 
[i. In the following, we fix the reference volume to be Popp's one. Lemma [20] can be used to 
generalize the results to the case of a generic //-divergence. 

With a divergence and a gradient at our disposal, we are ready to define the sub-Laplacian 
associated with the volume form /t. 

Definition 21. Let [i G Q n (M), f G C°°(M). The sub-Laplacian associated with fi is the 
second order differential operator 

A/:=div(V/), (48) 

This definition reduces to the Laplace-Beltrami operator when [i is the Riemannian volume. 
As a consequence of Eq. (j4"3|) and the Leibniz rule for the divergence div(/X) = X f + / div(X), 
we can find the expression of the sub-Laplacian in terms of any local frame X±, . . . , X^: 

k k 

div(V/) =^2^(Xi(f)Xi) = ^2x i (X i (f))+div(X i )X i (f) . (49) 

i=l i=l 

Then 

k 

A = Y,Xf + &v(X i )X i . (50) 

i=l 
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Remark 5. Observe that the second order term of A, namely the "sum of squares" in Eq. (|50p . 
does not depend on the choice of the volume. Indeed, only the first order terms depend on 
it through the divergence operator, which changes according to Lemma [52] upon a change of 
volume. 

Remark 6. If we apply Proposition [19] to the horizontal gradient Vg, we obtain 

/ fAg» = -[ (V/,V 5 ) M , V/, 5 GC °°(M). (51) 

JM JM 

Then A is symmetric and negative on C^°(M). It can be proved that it is also essentially 
self-adjoint (see [14]). 

Now we prove a useful formula for the divergence associated with Popp's volume. Analogous 
formulae for /i-divergences are easily obtained by an application of Lemma 1201 

Lemma 22. Let X\, . . . , X n be a local adapted frame. Let div be the divergence associated with 
Popp's volume. Then, for i = 1, . . . ,n 



divX- 



_. m n \ 

-^Tr^X^.))^^ J,. (52) 



Proof. Let v G Q 1 (M), and X, Y £ Vec(M). The Lie derivative obeys Leibniz rule: 

Cx{v{Y)) = {C x v){Y) + v{C x Y). (53) 
Then, if v\ . . . , v n is the dual frame associated with Xi, . . . , X n 

n 

£xy = -J24i» l > ( 54 ) 
i=i 

which is the "dual formulation" of Eq. ([2]). By Theorem [1] Popp's volume is 

V = 1 v 1 A ... A v n . (55) 

yrLjdetB,- 

Then, for i = 1, . . . , n, 



+ 1 (C x y A ... A / + ... + z/ 1 A ... A £ Xl ^ n ) • (56) 



'UjdetB, 

Eq. (|52p now follows from the definition of divergence, Eq. (|54|) and Eq. ()56p . □ 
Finally, Corollary [2] is a straightforward consequence of Lemma [22] and Eq. (|50p . □ 

5 Volume preserving transformations 

This section is devoted to the proof of Propositions and 
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5.1 Proof of Proposition [5] 

Let (j) G Iso(M) be a (local) isometry, and 1 < i < m. The differential induces a linear map 

4>* : 8)^9 -»• ®*XV( ff ) • (57) 
Moreover 0* preserves the flag Pc.C I*" 1 . Therefore, it induces a linear map 

^:VyV^^V^ q) /V^ y (58) 
The key to the proof of Proposition [5] is the following lemma. 
Lemma 23. 0* and 0* are isometries of inner product spaces. 

Proof. The proof for 0* is trivial. The proof for 0* is as follows. Remember that the inner 
product on is induced by the surjective maps 7Tj : &T> — > T> l /T>' l ~ 1 defined by Eq. (|14p . 

Namely let Y G Pj/P* -1 . Then 

ll^lloj/oj" 1 = min {ll- z 'll®o 9 s.t. 7Ti(Z) = 1"} . (59) 
As a consequence of the properties of the Lie brackets, Xj o <f>* = o Xj. Therefore 

ll^llo* iv^ 1 = ™ n {\\4 > * z \\®v !t(q) s.t. iii^Z) = 4>*Y} = WfaYWpi /v i-i . (60) 

ql 1 vyH > 4,{q)l <j>(q) 

By polarization, eft* is an isometry. □ 

Since gr q (V) = (B^/D'^/V 1 ^ 1 is an orthogonal direct sum, 0* : gr q (V) — > gr^^V) is also 
an isometry of inner product spaces. 

Finally, Popp's volume is the canonical volume of gr (P) when the latter is identified with 
T q M through any choice of a local adapted frame. Since is equal to under such an 
identification, and the latter is an isometry of inner product spaces, the result follows. □ 



5.2 Proof of Proposition [|J] 

Let \i be a volume form such that </>*// = \i for any isometry <f> G Iso(M). There exists / G 
C°°(M), / ^ such that V = ffi. It follows that, for any G Iso(M) 

/ M = p = ^ = (/o^)^=(/o^) /J , (61) 

where we used the Iso(M)-invariance of Popp's volume. Then also / is Iso(M)-invariant, namely 
(ft* f = f f° r an Y ^ Iso(M). By hypothesis, the action of Iso(M) is transitive, then / is 
constant. □ 
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